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. We consider dihadron fragmentation functions, describing the fragmentation of a parton in two 

' unpolarized hadrons, and in particular extended dihadron fragmentation functions, exphcitly de- 

I pendent on the invariant mass, Mh, of the hadron pair. We first rederive the known results on 

^ , Mh-integrated functions using Jet Calculus techniques, and then we present the evolution equa- 

' tions for extended dihadron fragmentation functions. Our results are relevant for the analysis of 

experimental measurements of two-particle-inclusive processes at different energies. 
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I. INTRODUCTION 
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I The fragmentation of partons into hadrons has been studied in detail in semi-inclusive processes with one hadron 

■ detected in the final state, such as e^e^ annihilation, Semi-Inclusive Deep-Inelastic Scattering (SIDIS) or hadron- 

■ hadron coUisions. Factorization theorems (see, e.g., P, H, [1]) allow to separate perturbatively calculable short-range 

: coefficient functions from long-distance nonperturbative fragmentation functions D'\~^^{z), describing the "decay" of 
the hard parton i into an observed hadron h with fractional energy z, provided a hard scale is available. This is the 
, case of the e+e" hX reaction, where the coefficient function is known at least to 0{as) [1, HI, with the running 
II ' strong coupling constant. The same fragmentation functions, D\~^'^, occur in the factorized formula for SIDIS at 
Qh! 0{as) [1, [![, combined with the specific process-dependent coefficient functions at the same accuracy. For hadronic 
1^ ' collisions, factorization is usually assumed, but has not been proven yet. 
'"7^ ' When considering semi-inclusive processes with two detected hadrons in the final state, e.g., e+e^ — > hih2X, a new 
J> . class of fragmentation functions, the so-called Dihadron Fragmentation Functions (DiFF), needs to be introduced 
to guarantee factorization of all coUinear singularities From this perspective, DiFF are analogous to fracture 
. functions in the space- like regime 0, [1] . The DiFF evolution equations have been recently reanalyzed in Ref. Q 

■ and the cross section for e+e" /11/12A' has been calculated to 0{as) in Ref. [l^. At O(a^), the production of two 
hadrons hi, /12 with fractional energies zi, Z2 and belonging to the same jet, is described by a DiFF, j)^^^^^^ (^^i, Z2), 
i.e., the fragmentation of a single parton i into the two hadrons. At 0{as), hadrons produced in the same jet could 
either come from the fragmentation of a single parton into two hadrons or by the fragmentation of two coUinear 
partons, i and j, into single hadrons. This implies that evolution equations for DiFF contain an inhomogcneous term 
of the form D'C^^ ® Di^''^ [l|. 

All these studies focused on DiFF as functions of the energy fractions zi and Z2, integrated over all the other 
kinematical variables of the produced hadron pair, including their invariant mass Mh- However, the largest amount 
of experimental information related to DiFF consists of invariant mass spectra of hadron pairs produced in e^e~ 
annihilation Semi-Inclusive Deep-Inelastic Scatteri ng ( SIDIS) ^'^'^ proton-proton collisions [13, 

[H, In this paper, using the techniques of Jet Calculus [EllOl we deduce the evolution equations for DiFF with 
an explicit dependence on Mh- In analogy with Ref. [2l|, we address them as extended Dihadron Fragmentation 
Functions (cxtDiFF). 



*Eleetronie address; |ceccopieri@fis. u nipr . it | 

tElectronic address: mareo.radiei@pv .lnfn.lt| 

Eleetronie address: , alessandro. bacchetta@desy.de | 



2 



DiFF turn out to have important applications in polarization studies, since they can act as spin-analyzers of the 
fragmenting quark [2^ . In particular, the transverse polarization st of the fragmenting quark can be related to 
the azimuthal orientation of the plane containing the two hadron momenta Pi and P2 , through the mixed product 
Pi X P2 - ST- The strength of this relation is described by the DiFF h^''^'^^'^^ [23]. In SIDIS with transversely polarized 
targets, this DiFF appears in combination with the transversity distribution function 23, [25l. [26j . thus providing a 
way to constrain this elusive partonic distribution (for a review on transversity, see Ref. j27|). The HERMES [28| and 
COMPASS [2^ collaborations have reported preliminary measurements of the induced spin asymmetry (at (Q^) « 2.5 
GeV^). In the meanwhile, the BELLE collaboration is planning to perform the extraction of the fragmentation 
function in e+e^ annihilation [s^], but at the higher scale y/s w 10 GeV [3l|. The invariant-mass dependence of 
this fragmentation function is unknown and nontrivial, as shown, e.g., by model calculations psi. [26l [s^. Therefore, 
the study of the evolution properties of extDiFF is also timely. 

The paper is organized as follows. In Scc.[TTl using Jet Calculus we recover the inhomogeneous evolution equations for 
DiFF derived in Ref. iJi),]. In Sec. lIIIl following the same lines we deduce the evolution equations for the corresponding 
extDiFF. Finally, in Sec. Ml we draw our conclusions. 



II. INTEGRATED DIHADRON FRAGMENTATION FUNCTIONS 



The cross section at order 0(as) for the e+e — > hX process, where a hadron has momentum and energy 
fraction z = 2Ph ■ q/Q^ with respect to the center-of-mass (cm) energy = q^, can be written formally as 

i 

where D\^^ describes the "decay" at leading twist of the hard parton i into the observed hadron ft,, the sum on i 
running over all possible partons species i = q,q,g. The process-dependent coefRcient function cr* can be calculated 
and regularized in perturbation theory and it is known at least at Ojcts) The D\~'^ at 0{as), i.e., which absorbs 

all the coUinear singularities, can be accurately parametrized 33, SJli except for the z — > portion of phase space. 

The generalization of Eq. ([T|) to the process e+e^ hih2X is not straightforward, if one wants to cover the whole 
phase space accessible to (/ii,/i2). The differential cross section, again at 0{as), has been recently calculated in 
Ref. [lO| and reads, with obvious notations, 

= y Dr""' ® or'"' +ya'<g> or'"'''' , (2) 

dzidz2 /l^ i i ^ i ' ^ ^ 

where the DiFF jj'-^^^^^ contains information on the fragmentation, at leading twist and 0{as), of the hard parton 
i directly into the observed hadron pair /ii,/i2. At order ©(a^), the first term of Eq. ([2]) would correspond to the 
back-to-back emission of a parton and an antiparton, eventually fragmenting in the hadrons hi and /i2 belonging to 
two well separated jets. The second term would apply instead to the case where the hadron pair is produced very close 
in phase space and it is detected inside the same jet while the other jet is inclusively summed over. However, at order 
0{as) a new kind of coUinear singularities arises in the partonic cross section; it corresponds to the configuration 
where each hadron is obtained from the fragmentation of a single parton, the two partons being almost coUinear, i.e., 
with a very small relative transverse momentum rx- These l/r^ singularities cannot be reabsorbed in each D\~^'^, 
because they do not correspond to the back-to-back configuration. Hence, they must be reabsorbed in £)*^''i''2^ ^j^^g 
making the two terms in Eq. ^ indistinguishable [l3| . 

As a consequence, after integrating over rr the DiFF must satisfy the following evolution equation [l^ 



,Dl--'^^'^Hzi,Z2,Q') 



27r J u{\ — u) \u / \1 — u ' ' 



(3) 



where here, and in the following, a sum over repeated parton indices is understood. The first term in the right- 
hand side represents the usual homogeneous evolution for the DiFF jj'^^^^^'^ jn complete analogy with the case of 
single-hadron fragmentation: the probability for the parton i to fragment into the hadrons /ii,/i2, is affected by the 
probability of emitting a parton j with momentum fraction u through the Altarelli-Parisi splitting vertex Pji{u) [35| . 
listed in Eqs. (|A.ip - (|A.4p of the appendix for convenience. The second term is a new inhomogeneous contribution 
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that corresponds to the probabihty for the parton i to spht in the two partons j and k with momentum fractions u 
and (1 — m), respectively, each one fragmenting in one of the two observed hadrons. The PJj, (u) are the AltareUi-Parisi 
splitting functions without virtual contributions [131, again listed in Eqs. (jA.6|) - (jA.9P of the appendix for convenience. 
From this point of view, the situation is similar to the DIS case in the target fragmentation region, since the DiFF 
can be conceived as the time- like version of the fracture functions in the space- like domain 0, Q . 

In the following, we will make use of Jet Calculus 0, [lO] and recover the evolution equation ([3]) within this 
formalism. We will consider the semi-inclusive production of two hadrons, hi and /12, belonging to the same jet and 
neglecting the emission of wide-angle hard partons (and related jets). Therefore, we will not perform a fixed-order 
calculation of the e+e^ —>■ hih2X cross section. Rather, we will consider a parton i with a large virtuality eQ^ 
(0 < e < 1), which fragments in two hadrons hi and /12 inside the same jet. The virtuality can be reconstructed 
from the invariant mass of the jet by a suitable jet-finding algorithm [36] . The phase-space structure of coUinear 
singularities singled out in fixed-order calculations can be translated in Jet Calculus as a degeneracy in all possible 
competing mechanisms, which could realize the desidered final state configuration [1, H^. Thus, the cross section is 
the sum of all production mechanisms, as in Eq. ([2]). 

We use as evolution scale, instead of eQ^. In Leading Logarithmic Approximation (LLA), this substitution 
induces only subleading corrections and thus is fully justified within this approximation. Moreover, it is convenient 
to replace this variable with the evolution variable 



Y = ^— In 

27r/3o 



also named the evolution imaginary time 20]. In Eq. Q, is the renormalization scale and Po — {llNc~ 2N f) / {12tt) 
is the one-loop /? function with Nc, Nf, the number of colors and flavors, respectively. In LLA, the running of is 
taken into account at one loop by 



1 



where Ag^^^ is the infrared scale. Hence, the differential evolution length is just 



Let us define the variable y as 



■In 



(7) 



with Qq and two arbitrary scales, and introduce the perturbative parton-to-parton time-like evolution function 
E^j{x,y), which expresses the probability of finding a parton j at the scale Qg with a momentum fraction x of the 
parent parton i at the scale Q^. The function Ej{x,y) can be shown to satisfy standard evolution equations [gI. [20j 

^E]{x,y)= ['^E^(^,y)p,,iu), (8) 
dy u \u y 

that can be iteratively solved by using the initial condition 

E}{x,y)\y=o^S,,6il~x) , (9) 

with Sji the Kronecker symbol. The Ej{x, y) resums leading logarithms of the type a" \n^{Q'^ /Qq), which show up in 
the collinear limit of perturbative calculations at the partonic level. As a crosscheck, we can expand Eq. ([S]) at order 
0{as) with the initial condition ([9]); neglecting for simplicity the running of as in Eq. ([7]), we get 



E]{x, y) = E]{x, Ql Q') « ^,,,5(1 + ^ P^x) In ^ . (10) 







Leading logarithmic contributions are therefore automatically accounted for at all orders through the function E. 
Consider now the fragmentation process i — > hih2X where hi and /i2 are detected within the same jet and with 
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relative transverse momentum < Q^. The corresponding cross-section, normalized to the jet-cross section, can be 
written as 



ajet dzidz2 

dw 



Z1+Z2 



2-1 \—.—,yo)E{w,Y-yo) 



yo Jz^ + Z2^ u{l-u) ^ \wu J \w{l-u) 



(11) 



The first term "A" is the convolution of the DiFF at some arbitrary (but still perturbative) factorization scale j/q with 
the parton-to-parton evolution function E'j . The second term "_B" , instead, represents the two separate single-hadron 
fragmentations, integrated over all possible generic intermediate scales y at which the branching j kl at partonic 
level might occur. Integration limits in Eq. (jlip are fixed by momentum conservation. Both terms are depicted in 
Fig.m In order to recover Eq.([3]), we now take the derivative of Eq. (fTTj) with respect to the variable Y. Using Eq. ([8]) 



1, Y, i 




FIG. 1: Double- and single-hadron fragmentations in Eq. The momentum fractions are indicated along with the scale 

and the parton indices. The black dot represents the parton evolution function E. 



and the definition of D^^''^''^ we get for the term 



±D\-:J:^^^{z,^Z2.Y) 



"'21+22 ^ \w w / dY 

f '^,d\-/^-^(^A.y)pm. 



(12) 



The derivative d/ dY of the "-B" contribution in Eq. (jlip produces two terms, since there is an explicit K-dependence 
in the upper integration limit. Using Eq. ^ and the same procedure as before, we get 

±Dr:-^-^i.,^z2^Y)^ f ^M(i--) ^^,PLiu)D\-^^ (^.y)d[--^ i^^V^ 

dY • Jzi+Z2^ u(l~u) \wu / \w[l — u) 



^dyf -^T" ^^,PL{u)±E)i^,Y^y)D\--^(^^^^ (13) 
/„ir„ w /^i Mil — u dY ■' \wu / \w{l — u] I 



yo J Z1+Z2 



Z2 



U{\ — U) \U / \\ — U J J ZI+Z2 ^ \U U 

Summing up Eqs. (fT2|) and (|13p. we get 

±Dr'^^'^^{z,,z2.Y) = ^i^^-'H'^. (fl,fa,y) p,,(,) 



dY Jzi+Z2 " \u u 



1-22 X (14) 



u{l — u) 
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which is exactly Eq. ([3]), after changing the scahng variable Y back to the more familiar Q^. This derivation is useful 
to adjust the formalism of Jet Calculus from Refs. [13] to the calculation of the e+e" — > hih2X cross section at 
0{as) from Ref. However, we want to stress again that even if the two expressions are formally identical, they 
have been derived from rather different approaches. Eq. ^ from Ref. ^ld\ applies to the full phase space for the 
production of two hadrons. Instead, Eq. (jl4[) gives the evolution of the parton i in a jet within which we identify the 
two detected hadrons hi and thus, it is valid only in the portion of phase space defined by the jet. We recall that 
if is the cm energy of the e~^e~ annihilation, the event is characterized by a certain number of jets resulting also 
from a large-angle hard parton emission; this is not included in the LLA used here. Therefore, the evolution scale in 
Eq. ([M]) must be intended as the invariant mass of the considered jet, i.e., cQ^ with < e < 1. 



III. EXTENDED DIHADRON FRAGMENTATION FUNCTIONS 



In the previous section, we considered the inclusive production of two hadrons hi and h2 inside the same jet, summing 
up all possible values of their invariant mass Mh- If the process starts from the hard scale of the fragmenting 
parton (or, equivalently, eQ^ in the case of the jet), there is no intermediate scale that allows to distinguish the two 
contributions in Eq. the j/q i^i the "A" term is arbitrary, and the scale y for the partonic branching in the "B" 

term is summed over. 

However, most of the experimental information on unpolarizcd DiFF consists of invariant mass spectra of hadron 
pairs [lllliirliljQEIiijEHIill- in addition, effects related to the partial- wave expansion of DiFF [13 are best 
explored when the latter explicitly depend on M^. Hence, in this section we will address the evolution equations for 
DiFF at the fixed scale M^. We will indicate these objects as extended Dihadron Fragmentation Functions (cxtDiFF), 
as the time- like analogue of the extended fracture functions that were introduced in Ref. [21] for the space- like SIDIS 
in the target fragmentation region. 

The dependence of the extDiFF upon can be easily mapped into R^, the square of the relative transverse 
momentum of the hadron pair. In fact, it results [37| 



{Pit - F2t)' 



Z1Z2 

Zl + Z2 



Ml 

Zl + Z2 



Ml 

Zl 



Ml 

Z2 



(15) 



with Ml and M2 the masses of the hadrons hi and ft,2j respectively. From the first line of Eq. pT|) . we get the obvious 
definition 



1 da 



i — >h\ h2 



a jet dzidz2 



i—^hi h'2 



D 



(zi, Z2, Y) = D\-X^^^izi,Z2,Y) + Dl-;^^^'^^ {zi,Z2,Y) 
1 da^^'^^^'^ f 2 



(Jjet dzidz2dR^ 



(16) 



J- = / dR'TD\-J'^"^{zi,Z2,R'T,Y)+ I dRlD\-;j^^''^{zi,Z2,R'T,Y) . 



The inhomogeneous "i?" term describes the time-like branching of parton i in two partons k and / with transverse 
relative momentum vt, eventually fragmenting in the two hadrons hi and h2 with transverse relative momentum Rt- 
If the R^ scale is fixed and in the perturbative regime, the scale at the partonic branching is no longer arbitrary as 
in Eq. ([TT|) [3, [l^]- In fact, if u and (1 — it) are the fractional momenta of partons k and I, the parton virtualities are 
related by 



1 



4u(l - u) 



(17) 



Hence, fixing rx at the partonic level determines in turn the branching scale kf. At the hadronic level, this is not 
guaranteed and some assumptions must be made. We will suppose that in the fragmentations k hi and / — > /12 
the parton virtualities are negligible, i.e., k^ ~ kf ~ 0, meaning that, once the branching i ^ kl has occurred, both 
the perturbative and nonperturbative transverse momenta generated in the fragmentation of the partons k and I are 
negligible (incidentally, perturbatively generated transverse momenta can be taken into account by using time-like 
evolution equations depending on transverse momentum [Hj). Consequently, the transverse relative momentum rx 
between k and / should not be substantially altered in the fragmentation, implying 



kf 



R. 



■T ■ 



(18) 



Corrections to the above relation affect our final result only at subleading level. Instead, the above assumption is 
also consistent with the approximation we are working with. Leading logarithms are known to manifest themselves 
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when the transversa momenta of the emitted partons are strongly ordered along the ladder, e.g., ^ ''t i ^ ■ • ■ ^ 
''t n Qg. Hadron pairs with large relative transverse momentum Rt are thus produced earlier in the imaginary 
time Y than hadron pairs with small Rt, as appropriate for time-like kinematics. 

If i?y is thus in the pcrturbativc domain, in analogy with Eq. ([4]), we can define the variable yx as 



Vt — In 



UsiR"^) 



or, in differential form, 

d as{R^) d 



(19) 



(20) 



27ri?| dyx ' 

Since the scale at which the branching occurs is fixed by i?^, from Eqs. (fTTj) . (fT6|) and ((20|) . we obtain 
D\-;J''''-{zi,Z2,R^T.Y) 

2ttR^ dyr Jyg Jz^+^^w^ Jzi u(l - u) \wu ) \w{\ ~ u) J (21) 



P,iiu) E] {w, Y - yr) dI^'^^ (^, y^) D[- 



2t:RI J,^+,^w^ u{1-u) ' ' ^ Kwu"') ^ V^"(l-") 



,yT 



If the scale i?^ is fixed in the nonperturbative regime, the above arguments leading to Eq. (|20p do not apply. This 
is the case for the homogeneous "A" term in Eq. pB]) . which describes the direct fragmentation of parton i in the 
two hadrons hi and h2'. the virtuality kf of the parent parton cannot be reconstructed from R'^ and it is set to the 
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arbitrary factorization scale Qq (or, in our notations, yo). From Eqs. (Ilip and p6p . we simply get 



Dr^^'^^{zi,z,,Rl,Y)= f ^Di'f'''^(-,-,R'T.yo) E]{w,Y- 



yo) . (22) 



Summing up Eqs. (|22p and (j2ip . and providing each term with an extra step function to separate the two different 
kincmatical regimes, we get the complete expression for the extDiFF at LLA: 

i?r''i''^(zi, Z2, Rl,Y)= D\-^/'''^zi, Z2, r) + D\J''''%zi,Z2,R^T, Y) 
dw 



D 

Z1+Z2 



^.-^/u/.. (^_1^'^^r2 y\ E](w,Y-y,)e{y,-yT) 



27ri?| J^^^^^ u;2 7ii u{l - u) 



a.iRl) dw du p.^^,^j,.^^^y_y^^j,.^..f^^^^ 

) ■' Vww / \w[l — u) J 

(23) 

Note that, despite the presence of the step functions, the separation between the two regimes is still arbitrary, since 
it depends on yo which is itself arbitrary. The evolution equations for the extDiFF can be obtained, in analogy with 
Eq. (fT4|) . by taking the derivative with respect to Y [or, equivalently, via Eq. ([6])]. By using Eq. ([8|) and the 
definition (1^5]) of the extDiFF themselves, we get 

dh^^i (zi,Z2,i?T,Q)-^;^y^^^^^^i^i [-,-,Rt,Q )PM^). (24) 

We explicitly checked that by integrating Eq. (IM)) upon R^ we recover Eq. . In conclusion, if the hadron pair is 
inclusively produced in the same jet at fixed transverse relative momentum Rt (or, equivalently, fixed invariant mass 
Mfi), the evolution equations for the extDiFF are of the standard homogeneous type. The explicit dependence on 
this new scale breaks the degeneracy of the two production mechanisms described in the previous section. The same 
arguments apply to the SIDIS target fragmentation region where extended fracture functions, explicitly depending 
upon the invariant momentum transfer t between the incoming and outgoing hadron, satisfy a homogeneous evolution 
equation [2lj. On the basis of Eq. ([M)) . we argue that the cross section at order 0{as) for the inclusive production of 
the two hadrons /ii, /i2, inside the same jet and with invariant mass Mh, can be expressed in the factorized form 

a'® ^^'.1/^2 (^2^^ g2) ^ (25) 



dzidz2dR^ 
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where cr' are the same coefficient functions found in the single-hadron inclusive cross section of Eq. ([T]). In our 
above derivation, we used the techniques of Jet Calculus 0, l2Clj |. where the factorization of coUincar singularities can 
be automatically accommodated through the use of the parton-to-parton evolution function E. Exchanges of soft 
particles are, however, not accounted for. 

Eq. ((24)) can be conveniently diagonalized using a double Mellin transformation. We define 

pi pl-zi 



and the anomalous dimension 



A) {n + m)= f du Pj, {u) u"+"-2 (37) 



With simple algebra manipulations, it is easy to verify that 



Kr,n'^' (Rt, Q') = Di7rt''' (Rt, Q') A] (n + m) . (28) 



dlnQ2 V >^ y 27r 

The above results can be extended also to^polarized extDiPF, in particular to the only one surviving when the 
adron pair is coUinear, i.i 
unpolarized case, namely 

d 



hadron pair is coUinear, i.e., jj^^^'^^'^^ [23|, |37[. The evolution equations for this function have the same form of the 



dhiQ: 



(29) 



where the splitting functions 5Pji for a transversely polarized fragmenting parton are used |39| . The 5Pji are listed 
in Eqs. (|A.10|) - (|A.13p of the appendix. 

These evolution equations can be conveniently used for phcnomcnological analyses, since they can connect experi- 
mental data taken at different energies. 



IV. CONCLUSIONS 



We have shown that in leading logarithm approximation the so-called extended Dihadron Fragmentation Functions 
(extDiFF), describing the inclusive production of two hadrons inside the same jet at fixed invariant mass M/i, satisfy 
evolution equations of the same homogeneous type as in the single-hadron fragmentation case. We stress that the 
explicit dependence on the scale is required to break the degeneracy at 0{as) between the fragmentation from 
a single parton or after the branching in two coUinear partons. While the first contribution pertains to the nonper- 
turbative regime, in the latter the transverse relative momentum Rt of the two hadrons can be traced back to the 
transverse relative momentum of the two coUinear partons after the branching, and, ultimately, to the hard scale of 
the originating parton. The analysis of the corresponding contribution to extDiFF shows that the dependence on this 
perturbative scale can be predicted. 

In our derivation, we used the techniques of Jet Calculus @, [20| . Factorization of coUinear singularities can be 
automatically accommodated through the use of the parton-to-parton evolution function E. On the basis of the 
simple result for the evolution equations of extDiFF, we argue that the cross section at order 0{as) for the inclusive 
production of the two hadrons /ii, /12, inside the same jet and with invariant mass M?i, can be expressed in a factorized 
form involving the same coefficient functions as in the single-hadron inclusive cross section. A complete proof of this 
statement would require however the inclusion of soft particle exchanges, which are not accounted for in Jet Calculus 
approach. 

Evolution equations can be extended also to polarized extDiFF, in particular to //^ '^''i^^ j2^^ Such fragmentation 
function can be extracted from e^e^ annihilation [30], evolved to the scale of semi-inclusive deep inelastic scattering 
measurements and allow the extraction of the transversity distribution function [23, [1^, [1^ ■ 
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APPENDIX: SPLITTING FUNCTIONS 



The unpolarized leading-order splitting functions P(u) read [35 



21^) +lsil-u)-l-u 



P,g{u) = Tr[1-2u + 2u'] 
'2 



PaM = Cf 
Pgg{u) = 2Ca 



-2 + u 



H h u(l - m) - 2 



IICa - 2N 



(5(1 -u), 



where Nf is the number of flavors, Ca = 3, Cp — 4/3, Tr = 1/2; the " + " prescription is defined as usual by 



dzg{z) [/(z)-/(l)]. 



(A.l) 

(A.2) 
(A.3) 

(A.4) 



(A.5) 



The unpolarized leading-order splitting functions P are readily obtained from the previous ones by dropping virtual 
contributions at the endpoint. In our notations, they read: 



1 - u 



pa I 



Tr{\-2u + 2u^\, 



Pg%{u) = 2Ca 



2 + u 

u 



1 1 / s ' 

H h w(l - u) - 2 



1 — u u 

The transversely polarized leading-order splitting functions 5P read (sol . |4i 



SPqgiu) = 0, 
SPggiu) = 0, 

= 2Ca 



1 - u 



Q -S{l-u). 



(A.6) 
(A.7) 
(A.8) 

(A.9) 



(A.IO) 

(A.ll) 
(A.12) 

(A.13) 



Due to angular momentum conservation, there is no mixing between quarks and gluons. The last splitting function 
can be used to evolve the DiFF for linearly polarized gluons, SG"^ [4l| . 
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